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Abstract 

A finite-difference time-domain (FDTD) numerical analysis is used to demonstrate that a toroidal 
solenoid, coaxial with an electric dipole, is a remarkable non-radiating configuration. It can be used 
to measure the dielectric permittivity of any ambient matter. It becomes a directional radiator 
at an interface between two dielectric media, depositing energy in the material with the highest 
polarizability. 
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I. INTRODUCTION 



Toroidal and supertoroidal structures are widely present in nature and a supertoroid 
was explicitly drawn by Leonardo in 1490. The simplest examples of such objects would 
be toroidal solenoids with currents in them. More generally, fractal complications of the 
simple toroidal wiring known as supertoroidal structures are discussed and toroidal arrange- 
ments of electric and magnetic dipoles have been discussed in the literature. Today the 
main biological journals feature an increasing number of papers on proteins, viruses and 
phages possessing elements of toroidal and supertoroidal symmetry. At the same time we 
witness a growing stream of theoretical papers on the electrodynamics and optics of toroidal 
and supertoroidal currents, toroidal nanostructures, toroidal microscopic moments and in- 
teractions of electromagnetic fields with toroidal configuration Q]. Recent studies of phase 
transitions in ferroelectric nanodisks and nanorods m and toroidal arrangements of spins 
in magnetic structures pj] show growing interest in studying toroidal structures from the 
materials research community. 

Here we report for the first time a rigorous finite-difference time-domain numerical anal- 
ysis proving that a toroidal solenoid with poloidal wiring coaxial with an electric dipole is a 
remarkable non-radiating configuration. The property not to radiate electromagnetic energy 
is based on the destructive interference between the fields created by each of its constituents. 
We show that this configuration may be used as a sensor for the dielectric permittivity of 
the ambient matter. It becomes a directional radiator at an interface between two dielectric 
media depositing energy in the material with highest polarizability. 

Non-radiating configurations are such oscillating charge-current distributions that do not 
produce electromagnetic fields in the radiation zone. An early work ^1 shows that the orbital 
motion of a uniformly charged spherical shell of radius R will not produce any radiation if the 
radius R of the shell is equal to lcT/2 where c is the speed of light, T is the period of the orbit 
and / is an integer number. The general problem for absence of radiation from an arbitrary 
localized charge distribution, exhibiting periodic motion with period T = 27t/ujo, has been 
addressed and it has been shown that such a system does not generate electromagnetic 
potentials in the radiation zone if the Fourier components J{luJor/cr, Iujq) are not present in 
the spectrum of the current density J{r,t). This criterion also explains the behavior of an 
orbiting uniformly charged sphere. It has been pointed out in [5] that this condition may not 
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be necessary. It indeed ensures the disappearance of the electromagnetic potentials in the 
radiation zone, however calculations of the power emitted by the system show that its value 
is zero provided that J{lujQr/cr,lujQ) oc r, which is a weaker sufficient condition. Indeed, 
the latter condition only requires the absence of the components transverse to the wave- 
vector. It has been proved rigorously 6] that the absence of the transverse components 
of the Fourier spectrum of the current density is a necessary and sufficient condition for 
absence of radiation. Interestingly, such a condition has appeared in an earlier study, 3|j 
in connection with electromagnetic self-force action and self-oscillations of a non-relativistic 
particle. 

The important conclusion that can be drawn from the earlier results is that two types of 
non-radiating configurations can exist in principle. For the first type the Fourier components 
J(ujr/cr, uj) of the current density are zero. Numerous examples of systems pertaining to this 
sort of non-radiating configurations exist - |^,|5|,|8|. A characteristic feature of these systems 
is that both the electromagnetic fields and the electromagnetic potentials are zero. For the 
second type of non-radiating configurations the Fourier spectrum is purely longitudinal i.e. 
J{ujr/cr, uj) oc r. Here the electromagnetic fields are zero but as we show the electromagnetic 
potentials may be finite. 

Interestingly, it is pointed out that the case of J{ujr /cr,uj) oc r corresponds to triv- 
ial spherically symmetric radial oscillations of the charge density. Nevertheless non-trivial 
examples can be created using toroidal structures. Recent papers and show that a non- 
radiating configuration can be constructed by combining an infinitesimal toroidal solenoid 
with poloidal current fiowing in its windings (i.e. along the meridians of the toroid) with an 
electric dipole placed in the center of the toroid. The explicit calculations of Q and P| show 
that while the electromagnetic fields disappear outside such a composite object, the electro- 
magnetic potentials survive. As we show here this particular structure belongs to the second 
type of non-radiating systems and that it is the longitudinal part of the Fourier-spectrum 
of the current density which is responsible for the residual electromagnetic potentials in the 
radiation zone. 

The results of Q and PI suggest that the non-radiating configurations involving toroidal 
solenoids may have a number of interesting electromagnetic properties. These properties 
however have never been studied in proper detail. This is the aim of the present study. 

The physical nature of the problem is extremely well-suited to numerical modelling using 
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FIG. 1: Toroidal surface (a) and its cross-section with the z — y plane (b). The triangles show the 
direction of the surface current jp . 

the FDTD method jl^ which will be used in our analysis. It gives the possibility to address 
the electromagnetic properties of this specific structure consisting of a toroid coupled to a 
dipole in full numerical detail. In addition, an assessment of what possible applications such 
structures might have is given. 

Exact compensation between the fields generated by a toroidal solenoid and an electrical 
dipole takes place for infinitesimal objects only. It therefore seems plausible that assessment 
of the extent to which the properties of the infinitesimal non-radiating configurations are 
preserved by finite-dimensional counterparts should precede possible experimental designs. 
It is important also to determine what is the behavior of these structures under certain 
(non-trivial) perturbations. 

II. INFINITESIMAL TOROIDAL SOLENOIDS AND NON-RADIATING CON- 
FIGURATIONS 



The electromagnetic properties of toroidal solenoids and toroidal helix antennas are stud- 
ied . detail . .efe.e„ce. flQQQSy Q. He. o... *e result. tKat wU. be used 

in our exposition are briefiy summarized. 

The current fiowing along the meridians of a toroidal solenoid (known also as poloidal 
current, Fig. 1) can be presented in the form (see e.g. isl): 



jp = V X M, 



(1) 



since V.jp = 0. In Eq. (H)) jp is the current density vector and M = (0, M^, 0) is the 
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azimuthal magnetization vector, given by 

NI(t) 

if (p — + < and zero otherwise. In Eq. Q is the total number of windings and 
I[A\ is the magnitude of the current. Pursuing this idea a step further the magnetization 
M can in turn be expressed as 

M = V X T, (3) 
where T = (0, 0, T^) is called toroidization vector. The general problem is difficult to perform 



analytically 
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14l | and any limit that preserves the correct properties while at the 
same time giving valuable mathematical simplification is a step worth taking. Such a step 
is c? ^ 0. This is a useful step because it gives the toroidization vector in the following form 
(see e.g. QQ). 

T. = 6'ir). (4) 

Assuming monochromatic time- dependence, oc exp{—iujt) , and using (4) the magnetic field 
created by the toroidal solenoid can be obtained in the form 

NIdB?k\ 1 

Hp = -[ik ){r X n)exp{ikr), (5j 

8 r 

where n is a vector of unit length pointing along the 2;-axis and k is the wave vector. 

A dipole can be introduced at the center of the toroid. If this dipole is modeled as a 
piece of wire of length L^i with the current strength being equal to Id the dipole moment 
amplitude, po, can be expressed through iLdId = ^^Po, where po = pon. The magnetic field 
of the dipole is Q] 



ujk . 1 , , .exp(ikr) 
Att ikr 



The time-averaged power P emitted by the composite object (an infinitesimal toroidal 
solenoid coupled to an electrical dipole) is given by 

where T = irNIdP? /2 and e is the relative dielectric permittivity of the ambient matter. 
This expression can be generalized to include higher-order multipole moments jl6| . 



Equation (|7j) can be rewritten in the form 
where 



(9) 



is the effective relative dielectric permittivity of the medium in which electromagnetic fields 
of the toroid and the electric dipole can compensate each other. This suggests that it 
should be possible to measure the relative dielectric permittivities of media (e.g. liquids) by 
adjusting experimentally the ratio of the currents 1^ and / until a minimum of the emitted 
power is detected. Then the relative dielectric constant of the material under investigation 
can be obtained from Q. 

It has been pointed out in that while the electromagnetic fields disappear when the 
compensation condition Q is satisfied the electromagnetic _potentials survive. However, 

nri 

there are examples of non-radiating configurations (see e.g. p, ISl]) for which both the elec- 
tromagnetic fields and the electromagnetic potentials are zero. The question is then what 
is the physical reason for that and what is the difference between both types of electromag- 
netic systems. Following 0, 0] it can be shown that the difference is in the current- density 
spectra. To see this consider the vector potential 



d-^r'. (10) 
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In the radiation zone the standard approximation 

Anr 

Now if the current density j{r,t) is expressed through its Fourier-transform 



18^ can be used and (fTO|l reduces to 
J jir',t - r/c + r.r' /cr)d^r' (11) 



Eq. (flTjl becomes 



j{r,t) = y"j(fc,cj)e-'("")c/=^fc(iu; (12) 



/j(^^,c.)e-(-^/^)da; (13) 



where Vq = r/r. As Eq. i\VS\\ shows only those components of the current density spectrum 
that correspond to |fc| = ui/c contribute to radiation ^]. It is immediately clear that if the 



condition 

l(^,c.)=0 (14) 

is satisfied then the vector potential vanishes. Using the continuity equation, the Fourier- 
components of the charge density can be expressed from the Fourier-components of the 
current density according to p{k,uj) = k.j{k,uj)/uj and by following a procedure similar 
to deriving Eq. (fT!?|) but this time for the scalar potential it can be shown that the scalar 
potential is also zero if (fT^ is satisfied. Therefore ()14p ensures that the electromagnetic 
system considered is a non-radiating configuration. This general statement is a sufficient 
condition P]. 

The results of j^l imply however that the condition (|TH) is not necessary. With the 
assumption of a monochromatic time-dependence (|T!?jl reduces to 

47rr c 

The electromagnetic fields can then be obtained using ff = V x A//io and E = iV x H/ueo. 
Beacause (jlSj) is valid in the radiation zone only, Vq and 1/r must be treated as constants 
in deriving the fields from the vector potential. The result is 

E = ^M'^^^r, xQx ro)e-''^(*-^/=) (16) 
y eo 47rcr 

and 

H = t'^^^^iro X j)e— (*-'^/-). (17) 

From (fT^ and ()17|) it is clear that the time-averaged Poynting vector, {S) = ^E x H*, can 
be presented in the form 

(5)oc|rox(Jxro)|Vo (18) 

The quantity vq x [j x r^) is the radiation pattern of the system. As can be seen from ()18jl 
a charge-current distribution will not emit electromagnetic energy if 

j± = roX (j X ro) = (19) 

which is a weaker sufficient condition compared to p4j) . The fact that (fTIHl is also a necessary 
condition for the absence of radiation can be seen by setting E and H to zero in ()16D and 
(fT7|) and this was e.ab.i..ed in 0. Howeve. « ^a. appealed in ..e eaHie. .udie. QQ. 



The identity j = vq x (j x Vq) + ro(r'o.j) = j_L + j\\ and the comparison of (fTI?|) with 
(jl4p show that systems satisfying (|14p will emit no electromagnetic energy and will not 
produce electromagnetic potentials. On the other hand for systems satisfying the weaker 
condition (fT^ the electromagnetic potentials are not necessarily zero, since the longitudinal 
(parallel to Vq) part of the current-density spectrum j\\ will contribute to the vector potential 
according to Eq. (fr5|) . 

It is easy to show that the non-radiating configuration consisting of a toroidal solenoid 
coaxial with an electric dipole is an example of the second type of non-radiating configura- 
tions - systems satisfying p9|) and not (fT^ . The current density associated with this system 
is j{r) = TV X (V X n6^{r)) + IdLdn6^{r) and its Fourier spectrum is given by 

(271)^ j{k) = -Tk{k.n) + {k^T + IdLd)n (20) 

For values of A; = uj^Tq/c and using the compensation condition (jH)), Eq. reduces to 

(2vr)^j(i^) = -T^n{r,.n). (21) 

c & 

As shows the current density spectrum is purely parallel to for wavenumber values 
|fc| = Lo^jc. It can be concluded that it is the survival of the longitudinal part of the current 
density spectrum that gives the possibihty to create non-zero electromagnetic potentials in 
the radiation zone in the absence of electromagnetic fields. 

III. NUMERICAL MODELING OF THE INTERACTION OF A NON- 
RADIATING CONFIGURATION WITH THE INTERFACE BETWEEN TWO MA- 
TERIALS 

As can be seen from Section 2 the composite emitter - toroid and dipole - becomes a 
non-radiating configuration (note that the compensation condition (jH)) is satisfied) due to 
the destructive interference between the electromagnetic fields created by the toroid and 
by the electric dipole. This interference occurs in all possible directions in a homogeneous 
medium. In an inhomogeneous material however as would be encountered for an interface 
between two dielectrics with relative permittivity constants ei, and 62 this does not happen. 
To appreciate this consider the situation shown in Figure El in which an emitter consisting of 
an electric dipole and a toroidal solenoid is placed in a medium with dielectric permittivity 
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PML 




FIG. 2: Non-radiating configuration consisting of a toroidal solenoid and an electric dipole near 
the interface between two materials. Cylindrical symmetry is assumed. PML - perfectly matched 
layer; d - dipole, T- toroid; Si, S2 - cylindrical surfaces used to calculate the power emitted in 
each of the materials with e = ei and e = €2, respectively. 



equal to ei. This medium is separated from a second one by a planar interface located 
at a distance D from the equatorial plane of the toroid. In the absence of the interface 
the system is non-radiative and the effective permittivity is e'= ei. In order to assess the 
consequences flowing from the presence of the interface between the two dielectrics and also 
the finite size of both the toroid and the dipole it is necessary to solve the Maxwell's equation 
exactly. This can be achieved numerically using the FDTD method 1^1 • The latter can be 
considerably simplified since the toroid is a body of revolution (BOR). Taking advantage of 
the axial symmetry reduces the problem to a two-dimensional one. Cylindrical coordinates 
can be used and there is no dependence on the azimuthal variable angle (Fig. 12)). This 
implementation of the FDTD method is known as BOR-FDTD |lC|. The computational 
domain is terminated by a standard perfectly matched layer (PML) 3, 3- The radiation 
of both the toroid and the dipole is categorized by the field components {Ep, E;^, H^) that 
are not identically zero and hence it is of E-type (TM) The applicability of the FDTD 
method to radiating structures (antenna problems) is well established and this technique 
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has been successfully applied to various designs j20|, |21|, |22|. In the model the poloidal 
current jp, is expressed through the azimuthal component of the magnetization which is 
consistent with the assumption that all the parts of the toroid respond simultaneously (or 
with negligible delay) to the driving voltage. This is expected to occur when the size of the 
toroid is much smaller than the wavelength. To evaluate the directional properties of the 
emitting structure studied, the quantities Pi and P2 are introduced and defined as 

Pi= [ (5).dSi,z = l,2 (22) 



In fj22j) (S) is the time-averaged Poynting vector and Si, S2 (see Fig. ^ are cylindrical 
surfaces placed away from the source (close to the PML region) in order to ensure that the 
near-field contributions have negligible effect on the power values calculated according to 

(1221). 



IV. RESULTS AND DISCUSSION 



In order to model the behavior of the non-radiating configuration the following parameter 
values are selected. The larger and the smaller radii of the toroidal solenoid are fixed to 
d = 1 cm and R = 0.5 cm, respectively, the dipole length is Ld=0.9 cm and the excitation 
frequency is uj/2tt = 1 GHz. The FDTD-grid resolution is Ap = Az = A/300, where A is 
the free-space wavelength. Since Eq. (jHl) is strictly valid for infinitesimal objects only, it is 
necessary to make sure that for the selected values of the parameters the contributions from 
the higher-order multipoles are negligible. To verify this Eq. (jSJ has been compared with 
results obtained from FDTD simulations in a homogeneous material (this pertains to the 
case of ei = €2 = e in Fig. |21) and the result is presented in Figure El The simulations are in 
good agreement with Eq. dHJ. This means that for the selected values of the parameters the 
contributions of the toroidal dipole moment and the electrical dipole moment are dominant. 

Figure 4 compares the directional properties of a perturbed non-radiating configuration 
with that of an electric dipole. The ratio between the power values Pi and P2 emitted 
in the materials with dielectric constants ei and €2, is computed for several values of the 
distance D using either a non-radiating configuration or an electric dipole. The presence of 
the interface affects both types of emitters. However Fig. HI shows that a larger fraction of 
the total emitted power can be directed in the material with e = 62 for the case in which the 
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FIG. 3: Normalized emitted power P„ = 12Pttc/ {fio{IdLiiio)'^) versus the relative dielectric per- 
mittivity e of the ambient dielectric material. The value of Z (Eq. ©) is ?= 1.5. The solid curve 
and the solid squares are the analytical result (Eq. Q) and the numerical result, respectively. 
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FIG. 4: The ratio between the powers Pi and P2 emitted in the materials with dielectric constant 
ei and €2, respectively, by a non-radiating configuration (NRG, solid squares) and an electric 
dipole (Electric dipole, solid triangles) as a function of the distance D between the emitter and the 
interface. The parameter ?of the non-radiating configuration is ?= ei. 

emitter is a non-radiating configuration. Comparing the performance of the non-radiating 
configuration with that of the dipole acting along shows that the non-radiating configurations 
has a clear advantage in the ability to direct a larger fraction of the total emitted power in 
a material with higher value of the dielectric constant. This advantage disappears in the 
proportion of the increase of the distance to the interface. It has been verified that the 
dependence of the ratio P2/P1 on the distance D for both the emitters stems mainly from 
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the dependence of the quantity P2 on D. The value of Pi appears to be less susceptible to 
the variations of D for this range of parameter values. For some applications it might be 
desirable to direct electromagnetic energy within a certain material while little or no energy 
is emitted to the surrounding space. It seems that a non-radiating configuration with e' = ei 
may be suitable for this purpose. Relatively far from the interface it does not radiate at all, 
or radiates a small amount of power. Bringing the non-radiating configuration into contact 
with the interface will lead to an increase of the total emitted power Pi + P2, (keeping the 
values of the currents Id and / fixed) with the contribution P2 predominating strongly. 

To study this property further Figure El and Figure El show the dependence of the ratio 
P2/P1 on the dielectric constant of the substrate for two fixed values of the distance D 
between the emitters and the interface and 'e = ei. The directional properties of the non- 
radiating composite object are compared with those of its constituents - the electrical dipole 
and the toroidal solenoid. As Fig. Eland Fig. Elshow, the ratio P2/P1 for the toroidal solenoid 
and for the electrical dipole shows little dependence on the dielectric constant of the substrate 
62. At same time, when non-radiating configuration is used as an emitter, not only the ratio 
P2/P1 is higher, but it increases strongly with the increase of €2- This shows that in the region 
of parameter values studied the directional properties of the non-radiating configuration 
improve with the increase of the contrast between the relative dielectric permittivities of the 
two materials. A comparison between Fig. El and Fig. El shows that as the non-radiating 
configuration approaches the interface its performance improves. Indeed, it can be concluded 
that the optimum performance is achieved when the non-radiating configuration is in direct 
contact with the interface. This feature is in agreement with Figure IH The dependence of 
the emission properties of the non-radiating configuration upon the values of the dielectric 
constant of the substrate is suitable for sensor applications. 

The results presented in Fig. El are visually presented in Fig. [71 and Fig. El where the 
spatial distribution of the time-averaged Poynting vector around two of the studied emitters 
- non-radiating configuration and a toroidal solenoid - is shown. The Poynting vector values 
are normalized to the value of the total emitted power to allow a comparison at identical 
total emitted powers to be made. As can be seen the electromagnetic field created by the 
non-radiating configuration in free space (e = 1) is weak compared to that emitted by the 
toroidal solenoid. This shows that a larger fraction of the total emitted power is deposited 
in the substrate material. 
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FIG. 5: The ratio between the powers Pi and P2 emitted in the materials with dielectric constant 
ei and 62, respectively, by a non-radiating configuration (NRC, solid squares), electrical dipole 
(Dipole, solid circles) and the toroidal solenoid (Toroid, solid triangles) as a function of the dielectric 
constant 62. The parameter e of the non-radiating configuration is ? = e = 1 and the distance 
between the emitters and the interface is D = 14.95 mm. 
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FIG. 6: The same as in Figure El but for D = 5.98 mm. 
V. CONCLUSIONS 



In conclusion we studied a remarkable non-radiating configuration consisting of a toroidal 
solenoid coupled to an electrical dipole. The property not to radiate electromagnetic energy 
is based on the destructive interference between the fields created by each of its constituents. 
We show that the interference effect depends on the dielectric characteristics of the ambient 
matter and the configuration may be used in dielectric permittivity measurements. It be- 
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FIG. 7: (Color online) Spatial distribution of the time-averaged and normalized Poynting vector 
modulus ( ^/ Sp + S"^) /{Pi + P2) for a non-radiating configuration. The distance between the 



emitter and the interface is D = 5.98 mm. Note that logarithmic scale is used for the values of the 
Poynting vector. The values of the other parameters are e = ei = 1 and €2 = 8. 
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FIG. 8: (Color online) The same as in Figure [3 but the emitter here is a toroid. 



comes a directional radiator at an interface between two dielectric media depositing energy 
in the material with the highest polarizability. 
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